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\D ■ ABSTRACT 
O 
O 

^ • In many magnetized, dilute astrophysical plasmas, thermal conduction occurs 

Qj . almost exclusively parallel to magnetic field lines. In this case, the usual stability 

^ . criterion for convective stability, the Schwarzschild criterion, which depends on 

^ . entropy gradients, is modified. In the magnetized long mean free path regime, 

instability occurs for small wavenumbers when {dP/dz){d\nT/dz) > 0, which 
^ . we refer to as the Balbus criterion. We refer to the convective-type instability 

^ . that results as the magnetothermal instability (MTI). We use the equations of 

, MHD with anisotropic electron heat conduction to numerically simulate the linear 

growth and nonlinear saturation of the MTI in plane-parallel atmospheres that 



fields can damp the growth of the MTI and affect the nonlinear regime. The 
instability saturates when the atmosphere becomes isothermal as the source of 



O . 

\0 . are unstable according to the Balbus criterion. The linear growth rates measured 

^ , from the simulations are in excellent agreement with the weak field dispersion 

'q_|! relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic 

O 

^. 

d '. free energy is exhausted. By maintaining a fixed temperature difference between 

the top and bottom boundaries of the simulation domain, sustained convective 
^ ! turbulence can be driven. MTI-stable layers introduced by isotropic conduction 

are used to prevent the formation of unresolved, thermal boundary layers. We 
find that the largest component of the time-averaged heat fiux is due to advective 
motions as opposed to the actual thermal conduction itself. Finally, we explore 
the implications of this instability for a variety of astrophysical systems, such as 
neutron stars, the hot intracluster medium of galaxy clusters, and the structure 
of radiatively inefficient accretion fiows. 

Subject headings: accretion, accretion disks — convection — hydrodynamics — 
instabilities — MHD — stars: neutron — turbulence 
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1. Introduction 

In many dilute, magnetized astrophysical plasmas, the electron mean free path between 
collisions can be many orders of magnitude larger than the ion gyroradius. In this regime, 
the equations of ideal magnetohydrodynamics (MHD) that describe the fluid plasma must 
be supplemented with anisotropic transport terms for energy and momentum due to the near 
free-streaming motions of particles along magnetic field hnes (Braginskii 1965). Thermal 
conduction is dominated by the electrons compared to ions by a factor of the square root of 
the mass ratio. In a non-rotating system, it is sufficient to neglect the ion viscosity (Balbus 
2004). A fully collisionless treatment may be done using more complex closures such as 
Hammett & Perkins (1990). 

The implications of anisotropic transport terms on the overall dynamics of dilute as- 
trophysical plasmas is only beginning to be explored (Balbus 2001; Quataert, Hammett, & 
Borland 2002; Sharma, Hammett, & Quataert 2003). One of the most remarkable results 
obtained thus far is that the convective stability criterion for a weakly magnetized dilute 
plasma in which anisotropic electron heat conduction occurs is drastically modified from 
the usual Schwarzschild criteria (Balbus 2000). In particular, stratified atmospheres are 
unstable if they contain a temperature (as opposed to entropy) profile which is decreasing 
upward. There are intriguing analogies between the stability properties of rotationally sup- 
ported flows (where a weak magnetic fleld changes the stability criterion from a gradient of 
specific entropy to a gradient of angular velocity), and the convective stability of stratified 
atmospheres (where a weak magnetic field changes the stability criterion from a gradient 
of entropy to a gradient of temperature). The former is a result of the magnetorotational 
instability (MRI; Balbus & Hawley 1998). The latter is a result of anisotropic heat con- 
duction. To emphasize the analogy, we will refer to this new form of convective instability 
as the magnetothermal instability (MTI). The MTI may have profound imphcations for the 
structure and dynamics of many astrophysical systems. 

In this paper, we use numerical methods to explore the nonlinear evolution and satu- 
ration of the MTI in two- dimensions. We adopt an arbitrary vertical profile for a stratified 
atmosphere in which the entropy increases upward (and therefore is stable according to the 
Schwarzschild criterion), but in which the temperature is decreasing upwards (and there- 
fore is unstable according to the Balbus criterion, [dP/ dz){d\n.T /dz) < 0). We confirm 
the linear growth rates predicted by Balbus (2000) for dynamically weak magnetic fields 
and numerically measure the growth rates for stronger fields. We find that in the nonlinear 
regime vigorous convective turbulence results in efficient heat transport. Full details are 
published in Parrish & Stone (2005), hereafter PS. 

These results may have implications for stratified atmospheres where anisotropic trans- 
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port may be present. The most exciting potential application is to the hot X-ray emitting 
gas in the intracluster medium of galaxy clusters (Peterson & Fabian 2005; Markevitch 
1998). The hot plasmas are magnetized and have very long mean free paths along the mag- 
netic field lines; thus, they are a prime candidate for this instability. For example, the Hydra 
A cluster has T ^ 4.5 keV and a density of n ~ 10~^ — 10~^cm^^ giving a mean free path 
that's almost one-tenth of the virial radius. Other applications are to the atmospheres of 
neutron stars with moderate magnetic fields and radiatively inefficient accretion flows. 



2. Physics of the MTI 

2.1. Equations of MHD and Linear Stability 

The physics of the MTI is described by the usual equations of ideal MHD with the 
addition of a heat flux, Q, and a vertical gravitational acceleration, g. 
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where the symbols have their usual meaning with E the total energy. The heat flux contains 
contributions from electron motions (which are constrained to move primarily along fleld 
lines) and isotropic transport, typically due to radiative processes. Thus, Q = Qc + Q,r, 
where 

Qc = -Xcbb • VT, (5) 

Qr = -XiiVT, (6) 

where xc is the Spitzer Coulombic conductivity (Spitzer 1962), b is a unit vector in the 
direction of the magnetic field, and xr is the coefficient of isotropic conductivity. 

Some progress can be made analytically in the linear regime of the instability. I introduce 
two useful quantities. 
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With WKB theory, one is able to obtain a dispersion relation. The details of this process 
are to be found in §4 of Balbus (2000). The most important result of the hnear analysis is 
the instability criterion, 

,2 2 x'rdPdlnT ^ 

k^^A-—,^ ^<0' 8 

PX oz oz 

where v\ = B^/4np is the Alfvcn speed. In the limit of infinitesimal wavenumber, the 
instability criterion shows that any atmosphere with the temperature and pressure gradients 
in the same direction is unconditionally unstable, i.e. 

dPdlnT ^ 

We refer to the instability criterion Eqn.[9] as the Balbus criterion. The instability criterion 
of the magnetorotational instability (Balbus & Hawley 1998) can be written 

kV. + i^>o. (10) 

where Q is the angular velocity. The similarity between the MRl and the MTl is self- 
evident. Strong magnetic fields arc capable of stabilizing short- wavelenth perturbations in 
both instabilities through magnetic tension. 



2.2. Computational Method 

We use the 3D MHD code ATHENA (Gardiner & Stone 2005) with the addition of an 
operator-split anisotropic thermal conduction module for our simulations. Our initial state is 
always a convectively stable state {dS/dz > 0) in hydrodynamic equilibrium. We implement 
two different boundary conditions for exploring the nonlinear regime. The first boundary 
condition is that of an adiabatic boundary condition at the upper and lower boundaries, i.e. a 
Neumann boundary condition on temperature. This situation is ideal for single-mode studies 
of linear growth rates. In the Neumann boundary condition the magnetic field is refiected at 
the upper and lower boundaries, consistent with the adiabatic condition on heat flow. The 
second boundary condition fixes the the temperature at the upper and lower boundaries of 
the atmosphere, i.e. a Dirichlet boundary condition on temperature. This setup is useful 
for driven simulations where wc wish to study the effects of turbulence. In this boundary 
condition, the magnetic field is again reflected at the boundaries of the box, but heat is 
permitted to flow across the boundary from the constant temperature ghost zones. With 
these choices of boundary conditions, the net magnetic flux penetrating the box is constant 
in time as there is zero Maxwell stress at the boundary. If the net magnetic flux penetrating 
the box is initially zero, then in two dimensions the magnitude of the magnetic fleld must 
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decay in time as a result of Cowling's anti-dynamo theorem, but otherwise the saturated 
state is not affected. 



3. Results 

3.1. Single-Mode Perturbation and Qualitative Understanding of the MTI 

By examining a single mode-perturbation to the background state, we can gain an 
intuitive understanding of the physical mechanism of the instability. We begin by perturbing 
a convectively-stable, but MTl-unstable atmosphere with a very weak sub-sonic and sub- 
Alfenic sinusoidal velocity perturbation in a box with adiabatic boundary conditions. The 
evolution of the magnetic field lines is shown in Figure 1 for several different times. Notice 
in the upper right plot that at x ~ 0.03 that a parcel of fluid (as traced by the frozen-in 
field line) has been displaced upward in the atmosphere. As this parcel of fluid comes to 
mechanical equilibirum with the background state, it adiabatically cools. The magnetic 
field line now is partially aligned with the background temperature state, thus thermally 
connecting this parcel of fluid with a hotter parcel deeper in the atmosphere. As a result, 
heat flows along the fleld, causing the higher parcel to become buoyant. This buoyant motion 
causes the field line to be more aligned with the background temperature, increasing the heat 
fiux, and generating a runaway instability. 

It is instructive to examine the behavior of the temperature profile of the atmosphere. 
Figure 2 shows the horizontally averaged temperature profile of this run at various times 
(normalized to the sound crossing time). The initial state is a linear temperature profile 
decreasing with height. As the instability progresses, the temperature profile becomes more 
and more isothermal. Saturation occurs, as one would expect, when the temperature profile 
is almost completely isothermal, since the source of free energy has been depleted. More 
details are to be found in §4 of PS. 

3.2. Linear Growth Rates 

By following §4.3 of Balbus (2000) and using the Fourier convention, exp(at -|- ikx), 
one can derive a weak-field dispersion relation for the MTI as 
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where N is the Brunt- Vaisala frequency, the natural frequency of adiabatic oscillations for 
an atmosphere. With the single-mode perturbation simulations we are able to measure the 
growth rate for a variety of situations. Figure 3 plots the nondimensionalized growth rate 
versus wavenumbcr for theory (solid line) and the measured values from simulations (crosses). 
As can be seen these are in very good agreement. 

There arc essentially two ways to suppress the growth of the MTI. First, strong magnetic 
fields can exert tension that limits the growth and saturation of the instability. Second, 
isotropic conduction, as would result from radiative transport, can effectively short-circuit 
the thermal driving along field lines necessary for this instability to occur. For more detailed 
analysis, we refer the reader to §3 of PS. 

3.3. Nonlinear Regime and Efficiency of Heat Transport 

In order to assess the efficiency of heat transport in the magnetothermal instability, one 
needs to examine multimodc simulations seeded with Gaussian white noise perturbations 
and conducting boundary conditions at the top and bottom of the domain. The simplest 
such set-up results in narrow, unresolved boundary layers at the upper and lower boundaries, 
thus, making the heat flux difficult to measure accurately. As an alternative, we utilize a 
more physically relevant simulation. This setup involves an atmosphere that is convectively 
stable throughout, but MTI unstable only in the central region. The surrounding regions 
are stabilized to the MTI through the addition of isotropic conducitivity. As a result, the 
central unstable region is well-resolved. Figure 4 shows the evolution of magnetic field lines 
as this instability progresses. The magnetic field lines shown essentially track the central 
unstable region; however, the third panel clearly shows a plume of fiuid that is penetrating 
into the stable layer as a result of convective overshoot. This phenomenon is well-known in 
the solar magnetoconvection literature (Tobias, et al. 2001; Brummell, Clune, & Toomre 
2002). In three dimensions this type of behavior greatly amplifies the magnetic field in a 
local magnetic dynamo. 

More quantitative measurements can be made by comparing the time- and horizontally- 
averaged vertical heat fiuxes and breaking it down into Coulombic, radiative, and isotropic 
components. Figure 5 shows these quantities plotted as a function of time. At the mid- 
plane, the oscillatory advective fiux is clearly dominant at any given instant in time; how- 
ever, averaged in time the advective heat fiux contributes roughly | of the total heat fiux. 
The Coulombic fiux, which is relatively constant in time, contributes the remaining i. 
To determine the heat conduction efficiency of this instability, we compare it to the ex- 
pected vertical heat fiux across the simulation domain for pure uniform isotropic conduc- 
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tivity, namely, Qq 3.33 x 10~^. The time-averaged heat conduction at the midplane is 
{Qtot,bo%) ~ 3.54 X 10~^, which indicates that the instabihty transports the entire apphed 
heat flux efliciently. 

4. Conclusions and Application 

The most important conclusion of this work is that atmospheres with dS/ dZ < are not 
necessarily stable to convection. In fact, dilute atmospheres with weak to moderate magnetic 
fields can be convectively unstable by the Balbus criterion resulting in an instability that 
we call the magnetothermal instability. We have verified using MHD simulations that the 
measured linear growth rates agree with analytic WKB theory as predicted by Balbus (2000). 
For adiabatic boundary conditions, we find the saturated state is an isothermal temperature 
profile, corresponding to the exhaustion of the free energy in the system. For a driven 
instability with conducting boundary conditions, we find that the MTI efficiently transports 
heat, primarily by advective motions of the plasma in the vigorous convection that results. 

The most promising application of this instability is to clusters of galaxies. Structure 
formation calculations assuming ACDM cosmologies predict monotonically decreasing tem- 
perature profiles of the intracluster gas (Loken, et al 2002). Observations of clusters with 
Chandra, such as Hydra A (DeCrandi & Modlendi 2002), however, indicate essentially fiat 
temperature profiles. The intracluster medium is dilute, magnetized, and has a mean free 
path that could be as high as one-tenth of the cluster virial radius. It may be that these 
temperature profiles are representative of the saturated state of the MTI. This possibility 
will be explored in future work. 
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Fig. 1. — Snapshots of the magnetic field fines for a single- mode perturbation with adiabatic 
boundary conditions at various times during the evolution of the instability, (upper left) 
Inititial condition; (upper right) Linear phase; (lower left) Non-linear phase; (lower right) 
Saturated state. 
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Fig. 2. — Vertical profile of the horizontally- averaged temperature profile in the single mode 
case at various times. The initial state is a monotonically decreasing temperature profile 
with respect to height, and the final state is isothermal. 
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Fig. 3. — Solutions of the MTI dispersion relation in the weak field hmit for an atmosphere 
with dlnT/dlnS — —3. The axes are normahzed to the local Brunt- Vaisala frequency, N. 
The crosses are growth rates measured from simulations. 
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Fig. 4. — Snapshots of the magnetic field in the run with stable layers, (far left) Early hnear 
phase; (middle left) early non-hnear phase, (middle right) The MTI drives penetrative 
convection into the stable layers, and at late times (far right) magnetic fiux is pumped into 
the stable layers. 
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Fig. 5. — Time evolution of the horizontally-averaged heat fiux at the midplane and 80% 
height of the simulation domain in the run with stable layers. The total heat fiux (thick solid 
line) is subdivided into Coulombic (thin solid line), radiative (dashed lined), and advective 
(dotted line) components. The instantaneous total heat fiux is dominated by advective 
motions. 



